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ABSTRACT 

This talk reviews recent progress in formulating the dynamics of Kaon Physics, 
within the framework of the 1/Ac-expansion in QCD. 



1 Introduction 

In the Standard Model, the electroweak interactions of hadrons at very low en- 
ergies can be described by an effective Lagrangian which only has as active de- 
grees of freedom the flavour SU{?>) octet of the low-lying pseudoscalar particles. 
The underlying theory, however, is the gauge theory SU (3)c x SU (2) ^ x U{1)y^ 
which has as dynamical degrees of freedom quarks and gauge fields. Going from 
these degrees of freedom at high energies to an effective description in terms of 
mesons at low energies is, in principle, a problem which should be understood 



in terms of the evolution of the renormalization group from short-distances to 
long-distances. Unfortunately, it is difficult to carry out explicitly this evolu- 
tion because at energies, typically of a few GeV, non-perturbative dynamics 
like spontaneous chiral symmetry breaking and color confinement sets in. 

The suggestion to keep the number of colours Nc in QCD as a free pa- 
rameter was first made by G. 't Hooft 0) as a possible way to approach the 
study of non-perturbative phenomena. Many interesting properties have been 
proved, which suggest that, indeed, the theory in this limit has the bulk of the 
non-perturbative properties of the full QCD. The spectrum of the theory in 
the large-iVc limit consists of an infinite number of narrow stable meson states 
which are flavour nonets . This spectrum looks a priori rather different to the 
real world. The vector and axial-vector spectral functions measured in e"*"e~ 
hadrons and in the hadronic T-decay show indeed a richer structure than just a 
sum of narrow states. There are, however, many instances where one is only in- 
terested in observables which are given by weighted integrals of some hadronic 
spectral functions. Typical examples of that are the coupling constants of the 
effective chiral Lagrangian of QCD at low energies, as well as the coupling 
constants of the effective chiral Lagrangian of the electroweak interactions of 
pseudoscalar particles in the Standard Model. Some of these couplings are the 
ones needed to understand ii'-Physics quantitatively. In these examples the 
hadronic world predicted by large-iVc QCD provides already a good approxi- 
mation to the real hadronic spectrum. It is in this sense that I shall show that 
large-iVc QCD is a very useful phenomenological approach for understanding 
non-perturbative QCD physics at low energies; X-Physics in particular. 

2 The Chiral Lagrangian at Low Energies 

Typical terms of the chiral Lagrangian at low energies are 

-^eff = 7^0 {D^UD^^U^) + iiotr (U^Fu^^UFD +■■■ 
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where J7 is a 3 x 3 unitary matrix in flavour space which collects the fields 
of the low-lying pseudoscalar particles, (the Goldstone fields of spontaneous 



chiral symmetry breaking.) The first line shows typical terms of the strong 
interactions in the presence of external currents It ; the second line shows typical 
terms which appear when photons, W's and Z's have been integrated out in 
the presence of the strong interactions. We show under the braces the typical 
physical processes to which each term contributes. Each term is modulated 
by a constant: Fq, Liq,... C.gg,... which encode the underlying dynamics 
responsible for the appearance of the corresponding effective term. 

There are two crucial observations concerning the relation of these low 
energy constants to the underlying theory, that I want to emphasize. 

• The low-energy constants of the Strong Lagrangian, like Fq and Lio, 
are the coefficients of the Taylor expansion of appropriate QCD Green's 
Functions. For example, with li^LniQ'^) the correlation function of a left- 
current with a right-current in the chiral limit, (where the light quark 
masses are neglected,) the Taylor expansion 

- Q^TIlb{Q^)W-^^ = - 4Lio Q^ + O (Q^) , (2) 
defines the constants Fq and Liq. 

• By contrast, the low-energy constants of the ElectroWeak Lagrangian, 
like e.g. C and gg, are integrals of appropriate QCD Green's Functions. 
For example, including the effect of weak neutral currents 0) , 

c-^.r fi - t^ct) (-0^n,^(Q^)) . (3) 
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Their evaluation requires a precise matching of the short- distance and 
the long-distance contributions of the underlying Green's functions. 

These observations are completely generic, independently of the l/A'c-expansion. 
The large-iVc approximation helps, however, because it restricts the analytic 
structure oi the Green's functions in general, and ^lr{Q'^) hi particular, to be 
a sum of poles only; e.g., in large- A'c QCD, 



where the sums are, in principle, extended to an infinite number of states. 



There are two types of important restrictions on Green's functions like 
^lr{Q^)- One type follows from the fact that, as already stated above, the 
Taylor expansion at low euclidean momenta must match the low energy con- 
stants of the strong chiral Lagrangian. The other type of constraints follows 
from the short- distance properties of the underlying Green's functions which 
can be evaluated using the operator product expansion (OPE) technology. In 
the large- A'c limit, this results in a series of algebraic sum rules which restrict 
the coupling constants and masses of the hadronic poles. 

2.1 The Minimal Hadronic Ansatz Approximation to Large- A^c QCD 

In most cases of interest, the Green's functions which govern the low-energy 
constants of the chiral Lagrangian are order parameters of spontaneous chi- 
ral symmetry breaking; i.e. they vanish, in the chiral limit, order by order 
in the perturbative vacuum of QCD. That implies that they have a power 
fall-out in at large-Q^, ( e.g., the function ULRiQ^) falls as 1/Q®.) 

That also implies that within a finite radius in the complex Q^-plane, these 
Green's functions in the large- A^c limit, only have a finite number of poles. The 
minimal hadronic ansatz (MHA) approximation is fixed by the minimum num- 
ber of poles required to satisfy the OPE constraints. In the case of the ^lr{Q^) 
function, the MHA approximation to the large- Ac expression in Eq. (Q) results 
in the simple function 
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Q^l^Lum = ^0 (Q2 + M^KQ^ + Ml) ■ 

Inserting this in Eq. (H) gives a prediction to the low-energy constant C which 
governs the electromagnetic tt"*" — tt'^ = Am^r mass difference, with the result 

Am^ = (4.9 ± 0.4) MeV , MHA to Large- A^c QCD , (6) 

to be compared with the experimental value 

Am^ = (4.5936 ± 0.0005) MeV , Particle Data Book 6) . (7) 

The shape of the function in Eq. (^), normalized to its value at = 0, is 
shown in Fig. 1 below, (the continuous red curve.) Also shown in the same plot 
is the experimental curve, (the green dotted curve) obtained from the ALEPH 
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Figure 1: MHA (red curve,) versus ALEPH (green dots) and other predictions 

collaboration data 0) ; as well as the different shapes predicted by other analytic 
approaches. Let us comment on them individually. 

• The suggestion to use large-TVc QCD combined with lowest order xPT 
loops, was first proposed by Bardeen, Buras and Gerard in a series of 
seminal papers 0. The same approach has been applied by the Dortmund 
group IliP, in particular to the evaluation of e'/e. In this approach the 
hadronic ansatz to the Green's functions consists of Goldstone poles only. 
The predicted hadronic shape, as shown by the BBG, HKPSB line (black 
dotted), is a constant. 

• The Trieste group evaluate the relevant Green's functions using the con- 
stituent chiral quark model (CxQM) lirr 113' 0) . They have obtained a 

^See refs. & S 0^ and references therein. 



long list of predictions 1121'', in particular e'/e. The model gives an edu- 
cated first guess of the low-Q^ behaviour of the Green's functions, as one 
can judge from the CxQM-curve (green dashed), but it fails to reproduce 
the short-distance QCD-behaviour. 

• The extended Nambu-Jona-Lasinio (ENJL) model ^ has a better low- 
energy behaviour, as the ENJL-curve (blue dot-dashed) shows, but it 
fails to reproduce the short-distance behaviour of the OPE in QCD. Ar- 
guments to do the matching to short-distance QCD have been forcefully 
elaborated in refs. , which also have made a lot of predictions; a large 
value for e'/e in particular. 

In view of the difficulties that these analytic approaches have in reproducing 
the shape of the simplest Green's function one can think of, it is difficult to 
attribute more than a qualitative significance to their "predictions" ; e'/e in 
particular, which requires the interplay of several other Green's functions much 
more complex than Ulb{Q^). 

3 Applications 

The large- A'^c approach that we propose in order to compute a specific coupling 
of the chiral electroweak Lagrangian consists of the following steps: 

1. Identify the underlying QCD Green's functions. 

2. Work out the short-distance behaviour and the long-distance behaviour 
of the relevant Green's functions. 

3. Make a large-Nc ansatz for the underlying Green's functions. 

We have tested this approach with the calculation of a few low-energy observ- 
ables: 

• The electroweak AmTr mass difference © which we have already discussed. 
For a review, see e.g. ref. where earlier references can be found. 



• The hadronic vacuum polarization contribution to the anomalous mag- 
netic moment of the muon , with the result 

a^lwp = (5.7±1.7) X 10-*, (8) 

to be compared with an average of recent phenomenological determina- 
tions 

(6.949 ± 0.064) X 10-^ (9) 

• The 7r° e+e" and rj f^'^fJ-^ decay rates , with the results shown 
in Table 1 below, where 

"(^ -'4^^^ (10) 



Table 1: Summary of branching ratios results 



Branching Ratio (10|) 


Large A^c Approach 


Experiment EP 


i?(7r" ^ e+e-) x 10** 


6.2 ±0.3 


6.28 ±0.55 


R(t]^ fi+fi-) X 10^ 


1.4 ±0.2 


1.47 ±0.20 


R(t]^ e+e-) X 10** 


1.15 ±0.05 


< 1.8 X 10^ 



These successful predictions have encouraged us to start a project of a system- 
atic analysis of non-leptonic i^-decays within this large-iVc approach. So far 
we have completed two calculations of if-matrix elements within this large-iVc 
approach, which we next discuss. 

3.1 The Sx-Factor of K° - K° Mixing 

The factor in question is conventionally defined by the matrix element of the 
four-quark operator Qas=2{x) = {sLYdL){sLltidL){x): 

(i?"|QAS=2(0)|i^") = "^flMlBM . (11) 

To lowest order in the chiral expansion the operator QA5=2(a;) bosonizes into 
a term of O(p^) 

Qas=2{x) ^ -^5A5=2(a^) [{D^U^W]^^ [{D,.U^)U]^^ , (12) 
^See e.g. Prades's talk at KAON2001 © and references therein. 



with (7as=2(a*) a low energy constant, to be determined, which is a function 
of the renormalization scale /j, of the Wilson coefficient C(/i) which modulates 
the operator Qas=2{x) in the four-quark effective Lagrangian. The coupling 
constant gAS=2(/^), which has to be evaluated in the same renormalization 
scheme as the Wilson coefficient C{^) has been calculated, is given by an 
integral EH^ conceptually similar to the one which determines the electroweak 
constant C in Eq. (^). The invariant Bk defined as 



When comparing this result to other determinations, specially in lattice QCD, 
it should be realized that the unfactorized contribution to this result is the one 
in the chiral limit, ft is possible, in principle, to calculate chiral corrections 
within the same large--/Vc approach, but this has not yet been done. 

The result in Eq. (O) is compatible with the old current algebra predic- 

|on|\ ' — ■ 

tion which, to lowest order in chiral perturbation theory, relates the Bk 
factor to the ti^tt^ decay rate. In fact, our calculation of the Bk factor 

can be viewed as a successful prediction of the tt+tt" decay rate! 

As discussed in El) the bosonization of the four-quark operator Q\s=2 
and the bosonization of the operator Q2~Qi which generates A/ —1/2 transi- 
tions are related to each other in the combined chiral limit and next-to-leading 
order in the 1/Nc- expansion. Lowering the value of Bk from the leading large- 
Nc prediction Bk = 3/4 to the result in Eq. is correlated with an increase 
of the coupling constant gs in the lowest order effective chiral Lagrangian, (see 
Eq. (^,) which generates A/ = 1/2 transitions, and provides a first step to- 
wards a quantitative understanding of the dynamical origin of the A/ = 1/2 
rule. 




(13) 



can then be evaluated, with no free parameters, with the result EIP 



Bk = 0.38 ±0.11. 



(14) 



3.2 Electro Weak Four-Quark Operators 



These are the four-quark operators generated by the so called electroweak 
Penguin like diagrams ^ 

•••C7(m)Q7 + C8(m)08, (15) 

with 

Qt = Q{sLl''dL) V eq{qRjf,qB) and Qs = -12 V eg{sLqR){qRdL) , 



X! ^q^qRli^qR) and Qs = -12 ^ 

q—u,d,s q—u,d,s 



(16) 

and Cy (/i) , Cg (a*) their corresponding Wilson coefRcients. They generate a term 
of ©(p") in the effective chiral Lagrangian EH); therefore, the matrix elements 
of these operators, although suppressed by an factor, are chirally enhanced. 
Furthermore, the Wilson coefhcient Cs has a large imaginary part, which makes 
the matrix elements of the Qg operator to be particularly important in the 
evaluation oi e' /e. 

Within the large-A'c framework, the bosonization of these operators pro- 
duce matrix elements with the following counting 

{Q7)\o(pO) ^ 0{N,) +0{N",) and {Qs}\oipO) ^ 0{N^) + 0{N^ (17) 

Zweig suppressed 

A first estimate of the underlined contributions was made in ref. El) . The inclu- 
sion of final state interaction effects based on the leading large- A'c determina- 
tion of {Qs), (and (Qe),) in connection with a phenomenological determination 
of e'/e, has been recently discussed in E^) . 

The bosonization of the Qj operator to 0{p^) in the chiral expansion 
and to 0{Nc) is very similar to the calculation of the Z-contribution to the 
coupling constant C in Eq. (H). An evaluation which also takes into account 
the renormalization scheme dependence has been recently made in HI) . 

The contribution of 0{N^) to {Qs)\o{p°) is Zweig suppressed. It involves 
the sector of scalar (pseudoscalar) Green's functions where it is hinted from 
various phenomenological sources that the restriction to just the leading large- 
Nc contribution may not always be a good approximation. Fortunately, as first 



*See e.g., Buras lectures 



pointed out in , independently of large-iVc considerations, the bosonization 
of the Qs operator to 0{p^) in the chiral expansion can be related to the 
four-quark condensate (O2) = {0\{sLSR){dRdL)\0) by current algebra Ward 
identities, the same four-quark condensate which also appears in the OPE of 
the llLji{Q^) function discussed above. The crucial observation, here, is that 
large- A^c QCD gives a rather good description of the nL/j;((5^)-function, as we 
have seen earlier; in particular it implies that 

hm {-Q'llLB{Q'))Q' = J2fvMt^-T.fAMA- (18) 

•^^"^ V A 

This results in a determination of the matrix elements of {Qs)\o(p°) which does 
not require the separate knowledge of the Zweig suppressed 0{N^) term in 
Eq. (|l|). 

The numerical results we get for the matrix elements 

M7,8^((7r7r),.2|Q7,8l^") 2GeV (19) 

at the renormalization scale jj, — 2 GeV in the two schemes NDR and HV and 
in units of GeV'^ are shown in Table 2 below, (the first line.) 



Table 2: Matrix elements results, (see Eq. ^191)) 



METHOD M7(NDR) Af7(HV) MgCNDR) Af8(HV) 

Large— A^e Approach 

Rcf. H) 0.11 ±0.03 0.67± 0.20 3.5 ±1.1 3.5± 1.1 

Lattice QCD 

Ref. i^) 0.11 ±0.04 0.18 ±0.06 0.51 ± 0.10 0.62 ±0.12 

Dispersive Approach 

Ref. 0) 0.22 ±0.05 1.3 ±0.3 

Ref. H) 0.35 ±0.10 2.7 ±0.6 

Ref. H) 0.18 ± 0.12 0.50 ±0.06 2.13 ±0.85 2.44 ±0.86 

Ref. H) 0.16 ±0.10 0.49 ±0.07 2.22 ± 0.67 2.46 ± 0.70 



Also shown in the same table are other evaluations of matrix elements with 
which we can compare scheme dependences explicitly ^. Several remarks are 

There is a new "dispersive determination" in the literature 0) since the 
KAON2001 conference, but it is controversial as yet; this is why we do not 
include it in the Table. 



in order 



• Our evaluations of Ms do not include the terms of 0{a^) because, as 
pointed out in Ref. ^3) , the available results in the literature H) were 
not calculated in the right basis of four-quark operators needed here. 

• We find that our results for My are in very good agreement with the 
lattice results in the NDR scheme, but not in the HV scheme. This 
disagreement is, very likely, correlated with the strong discrepancy we 
have with the lattice result for M8(NDR). 

• The recent revised dispersive approach results , which now include 
the effect of higher terms in the OPE, are in agreement, within errors, 
with the large-TVc approach results. In fact, the agreement improves 
further if the new 0{a1) corrections, which have now been calculated in 
the right basis HI) , are also incorporated in the large-iVc approach. 

• Both the revised dispersive approach results and the large-iVc approach 
results for Mg are higher than the lattice results. The discrepancy may 
originate in the fact that, for reasonable values of {"tpip), most of the con- 
tribution to Mg appears to come from an OZI-violating Green's function 
which is something inaccessible in the quenched approximation at which 
the lattice results, so far, have been obtained. 
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Discussion 



Martinelli - Don't you think that your results for the Mg matrix elements 
could be affected by chiral corrections? 

de Rafael - There are chiral corrections, of course, which so far have 
not been evaluated; but the results that you reported in your talk, showing 
that the lowest order chiral Ward identities for these matrix elements are well 
reproduced in lattice QCD, are an indication that the chiral corrections should 
be rather small in this case. 

Isidori - Does your calculation of the P — s- 7*7* amplitude, used to predict 
P ^ II widths, reproduce well the available data on the P 7*7 form factor? 

de Rafael - Yes, but I suggest you to discuss this issue with Andreas 
Nyffeler, who has recently performed a detailed study of this subject with 
Marc Knecht in |hep-ph/0106034 . 

Roberts - Could your large- A^c methods be applied to a calculation of the 
hadronic light-by-light contribution to 5^ — 2? 

de Rafael - Yes, and in fact we have been discussing quite a lot about this 
in Marseille. We think that it should be possible to do a rather clean evaluation 
of the contribution which comes from two (PVV) three-point functions. 

Buras - I am puzzled by the fact that your hadronic integrals go all the 
way to infinity. 

de Rafael - This is due to the fact that we both are using a renormaliza- 
tion scheme based on dimensional regularization: you at short-distances and 
we at long-distances. We both depend on a scale /x which has to cancel in the 
matching. If we were doing a cut-off renormalization, the long-distance inte- 
grals would then be cut at a certain scale A indeed; but then you would have to 
redo all the short-distance analyses in the same cut-off renormalization which, 
as Donoghue pointed out, would introduce operators of higher dimension at 
short-distances. 



